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ABSTRACT 

The broad (FWHM - 10,000 km s"!) double-peaked Ha profile from the LINER/Seyfert 1 nucleus 
of NGC 1097 was discovered in 1991, and monitored for the following 11 years. The profile showed 
variations attributed to the rotation of gas in a non-axisymmetric Keplerian accretion disk, ionized 
by a varying radiatively inefficient accretion flow (RIAF) located in the inner parts of the disk. We 
present and model 11 new spectroscopic observations of the double-peaked profile taken between 2010 
March and 2011 March. This series of observations was motivated by the finding that in 2010 March 
the flux in the double-peaked line was again strong, becoming, in 2010 December, even stronger than 
in the observations of a decade ago. We also discovered shorter timescale variations than in the 
previous observations: (1) the flrst, of ^7 days, is interpreted as due to reverberation" of the variation 
of the ionizing source luminosity, and the timescale of 7 days as the light crossing time between the 
source and the accretion disk; this new timescale and its interpretation provides a distance between 
the emitting gas and the supermassive black hole and as such introduces a new constraint on its 
mass; (2) the second, of « 5 months, was attributed to the rotation of a spiral arm in the disk, which 
was found to occur on the dynamical timescale. We use two accretion disk models to fit theoretical 
profiles to the new data, both having non-axisymmetric emissivities produced by the presence of an 
one-armed spiral. Our modeling constrains the rotation period for the spiral to be w 18 months. 
This work supports our previous conclusion that the broad double-peaked Balmer emission lines in 
NGC 1097, and probably also in other low-luminosity active nuclei, originate from an accretion disk 
ionized by a central RIAF. 

Subject headings: accretion, accretion disks - galaxies: individual (NGC 1097) ~ galaxies: nuclei — 
galaxies: Seyfert - line: profiles 



1. INTRODUCTION 

The barred spiral galaxy NGC 1097 has an ac- 
tive nucleus classified as LI NER (Low - Ioniza tion Nu- 
clear Emission-Line Region, iHeckmanI ()1980l )). whose 
spectrum was found in 1991 to show transient 
broad (FWHM Rj 10,000 km s' ^) double-peaked Balmer 
lines (jStorchi-Bergmann et al.l [1993). The nucleus of 
NGC 1097 was the first low-luminosity active galactic nu- 
cleus (LLAGN) in a LINER discovered displaying such 
lines, and, because of these broad lines, has also been 
classified as a Seyfert 1 nucleus. Later, other LLAGNs 
in LINER nuclei were found to display double -peaked 
Balme r fines, such as M 81 (|Bower et al l I1996D. NGC 
4203 ([Shields et al.l [2000l) . and NGC 4450 (|Ho etall 
[200^ . 

The double-peaked Ha profile of the NGC 1097 nu- 
cleus was monitored for the following 11 years after the 
discovery (Storchi-Bergmann et al.l[2003f ) (SB03), show- 
ing variability in its flux, width and relative intensity 
of the blue and red peaks. These variations were at- 
tributed to gas emission from a thin Keplerian accre- 
tion disk, with a non-axisymmetric perturbation around 
a supermassive black hole (hereafter SM BH), where the 
gas is ionized by a central sou rce, (Eracl eous et al.lfl995l : 
iStorchi-Bergmann et al.lll995l[T997..200C 

The variability of the double-peaked profile of 



NGC 1097 is not rare. The first active nucleus found 
to display variable double-peaked broad Balmer emis- 
sion lines was the nucleus of the radio galaxy 3C 390.3 
(|Yee fc Okelll98lD . This object can be considered pro- 
totypical of the class of double-peaked emitters, because 
many of the now well-known features of these objects, 
such as variabi lity in the relative intensity of the blu e 
and red peaks (jVeilleux et al.l 119911 : IZheng et al.lll991| ). 
were first observed in 3C 390.3. Another remarkably well 
studied double-peaked emitter is the active nucleus of 
Arp 10 2B. For this object, a study bv iChen fc Halperiil 
(|1989f ) has shown that the double-peaked Balmer lines 
cannot be driven by local viscous dissipation in the line- 
emitting part of the disk, as the Ha luminosity is of 
the order or even exceeds the energy locally available to 
power the line. The same conclu sion was obtained for 
other doub le-peaked emit ters by 'E racleous fc HalpernI 
(|1994D and iStrateva et all (2006, 2003). These authors 
suggested that the lines are powered by illumination from 
an external source. The ionizing source could be an geo- 
metrically thick accretion flow, located within the inner 
radius of the line- emitting portion of t he disk. In the 
case of NGC 1097. lNemmen et all (|2006[ ) found that the 
spectral energy distribution (SED) of the nucleus is in- 
deed well described by such a st ructure, a radiatively in- 
efficient accretion flow (RIAF, iNaravan fc McClintock.1 
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(|2008D ). This has led to the conclusion that an in- 
ner RIAF is the ionizing source of the accretion disk in 
N GC 1097. 

iLewis et al] (|2010( ) and iGezari et aH (|2007D reported 
a long-term monitoring of 14 other double-peaked emit- 
ters and found similar variations to those observed for 
NGC 1097, 3C390.3 and Arp 102B. These objects also 
display in at least one epoch the red peak more in- 
tense than the blue peak. This indicates that most line- 
emitting disks are non-axisymmetric, as in an axisym- 
metric disk the blue peak should be always more intense 
than the red because of Doppler boosting. The disks 
could then be elliptical (jEracleous et al.|[l995D , or have a 
non-axisymm etric emitting stru cture, such as a hot spot 
or spiral arm ([Lewis et al.ll20ldl) . 

Another model which can reproduce double- 
pea ked emission-lin e profiles is the one proposed 
bv iGaskell. Khmek k Nazaroval (|2007[ ) and Gaskell 
(|2010[) . in which the ionizing source is the continuum 
from a thin accretion disk, and the origin of the double- 
peaked lines is the broad-line region (BLR) which is 
proposed to be an inner extension of the torus, in the 
shape of a "bird's nest". This geometry is not that 
distinct from that of the outer parts of an accretion disk, 
and in this regard is in approximate agreement with our 
proposed geometry. But , in th e case of NGC 1097, the 
study by iNemmen et al.l (|2006[ ) favors an RIAF origin 
for the nuclear ionizing continuum instead of a hot and 
thin accretion disk. 

The monitoring of double-peaked profiles is a valu- 
able experiment in the investigation of the structure 
and dynamics of the accretion disks and how they 
evolve as the gas moves inwar ds to feed the SMBH. In 
iStorchi-Bergmann et al] (|2003[ ). observing the NGC 1097 
Ha double-peaked profile approximately once or twice a 
year, we concluded that the asymmetry of the double- 
peaked profile was due to a spiral arm perturbation in 
the disk, rotating with a period of 5.5 years. But as rele- 
vant timescales for evolution of the innermost regions of 
AGNs can be as short as weeks and months, we decided 
to repeat our experiment on NGC 1097 in order to probe 
shorter timescales. Here we report the results of this ex- 
periment: we indeed found short variability timescales. 
We also verified that, although our last observations of 
2001 showed that the profile was becoming fainter and 
fainter, it is now again strong and was even stronger in 
2010 December than in all previous observations. 

This paper is organized as follows: in Section 2 we de- 
scribe the observations and the data reduction in Section 
3 we present the observational results, in Section 4 we de- 
scribe the models and the results from the modeling. In 
Section 5 we discuss the timescales of the accretion disk 
variability and the interpretation of the modeling. The 
conclusions of this work are presented in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

Long-slit optical spectra of the nucleus of NGC 1097 
were obtained in queue mode in 11 epochs from 2010 
March 4 to 2011 March 21 using the Gemini South Muhi- 
Object Spectrograph (GMOS) in the long-slit mode. The 
allocation was part of the "Poor Weather" opportunity, 
which is ideal for spectral monitoring which does not 
need absolute flux calibration or regular cadence. This 
is the case for the nucleus of NGC 1097, for which we 
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Observation Log 








p.a. 


Seeing 


Program ID 


UT Date 


(deg) 


(arcsec) 


GS-2009B-Q-99 


2010 Mar 4 


70 


1.00 


GS-2010A-Q-81 


2010 Aug 15 


55 


0.98 


GS-2010B-Q-90 


2010 Dec 22 


280 


0.63 


GS-2010B-Q-90 


2011 Jan 11 


280 


0.73 


GS-2010B-Q-90 


2011 Jan 15 


280 


0.73 


GS-2010B-Q-90 


2011 Jan 18 


280 


0.70 


GS-2010B-Q-90 


2011 Feb 8 


280 


1.21 


GS-2010B-Q-90 


2011 Feb 10 


280 


0.79 


GS-2010B-Q-90 


2011 Feb 14 


280 


0.72 


GS-2010B-Q-90 


2011 Mar 17 


280 


1.43 


GS-2010B-Q-90 


2011 Mar 21 


280 


0.82 



use the integrated flux in the narrow emission lines to 
normalize the spectra to match previous observations, 
under the assumption that the narrow lines do not vary 
within the extraction aperture. 

The first three observations were obtained at time in- 
tervals of 4-5 months, since our previous monitoring sug- 
gested significant variations on a timescale of years. How- 
ever, our new observations showed strong variations in 
the nuclear spectrum on timescales of months, which mo- 
tivated us to ask for new observations at shorter time 
intervals. We were then able to get at least two obser- 
vations a month from 2011 January to 2011 March, as 
shown in the observing log presented in Table [TJ This 
table lists the identification of the program, date of ob- 
servation, position angle (p.a.) of the slit and the average 
seeing during the observations, measured as the FWHM 
of a spatial profile extracted along a range of columns 
containing only the broad line (after subtraction of the 
underlying continuum) . 

The data were obtained with the GMOS South with a 
slit width of I'.'O {^ 80 pc), using the B600 grating and 
the GG455 filter to block possible higher order contami- 
nation. The spectral range covered was 6000-7000A for 
most spectra, at a spectral resolution of ~ 4.5A (FWHM 
of the lines in the arc spectrum, w 200 km s~^). Six 600s 
exposures were obtained at each epoch, which were then 
combined into a final spectrum. 

The slit width of I'.'O was selected because it is usu- 
ally larger than the average seeing during the observa- 
tions and was also the slit width of previous observations 
(SB03). Using the same slit width allows us to normal- 
ize the new spectra so that they match the flux of the 
narrow emission lines observed in previous spectra and 
study the evolution of the flux and profile of the broad 
Ha line over « 20 yr. 

The data reduction was performed using the standard 
tasks from the gmos package in the Image Reduction and 
Analysis Facility software (IRAFQ) . We extracted the nu- 
clear spectra using a window of I'.'O x I'.'O centered at 
the peak of the continuum emission, assumed to corre- 
spond to the galaxy nucleus. This position coincides also 
with the location of the unresolved source of the double- 
peaked Ha lines. Figure [1] shows the acquisition image 
of the observation of 2010 August 15, where the nuclear 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — Acquisition image from 2010 August 15, where we il- 
lustrate the locations of the I'.'Oxl'/O extraction windows for the 
nucleus and stellar population spectra. The square labeled A rep- 
resents the nuclear extraction window, while squares B and C rep- 
resent the extranuclear ones, from which we obtained the spectrum 
of the underlying stellar population. 

extraction region is represented by the square labeled A. 

2.1. Stellar population contribution 

Since the nuclear spectra show strong absorption lines 
from the underlying stellar population, we subtracted 
the stellar population contribution, in order to isolate 
and better study the emission line profile and in partic- 
ular the broad Ha line. A template stellar population 
spectrum was obtained as follows. For each epoch, in 
addition to the nuclear spectrum, we extracted two more 
spectra using two similar windows of I'.'O x I'.'O, one cen- 
tered at I'.'S to the southwest (Figure [H region B) and 
another at I'.'S to the northeast of the nucleus (Figure [H 
region C). These extranuclear spectra do not show any 
broad Ha emission, but show some narrow-line emission 
which was edited out, using as guid es, stellar popula - 
tion templates from previous studies (jBica et al.l 119881 ) . 
Averaging these two spectra, we obtain a stellar pop- 
ulation spectrum from a location approximately 100 pc 
away from the nucleus. Under the assumption that this 
spectrum is representative of the stellar population at 
the nucleus, we subtracted it from the nuclear spectrum 
(after normalizing the continuum to the corresponding 
value at the nucleus) thus isolating the gas emission. 

Figure[2]illustrates the process for the spectrum of 2010 
August 15, which covers a broader spectral range than 
most spectra, from 4750 to 7200A. We note that this 
process reveals the presence of a broad double-peaked H/3 
line which was not obvious prior to the subtraction, as 
well as some weak emission lines such as [N i] A5199A and 
[O i] A6363A. The Na I D absorption line in the residual 
spectrum is most probably interstellar, as all other stellar 
absorption lines did not show any residual. 

We also note that after subtracting the stellar popula- 
tion spectrum from the nuclear one, we do not see any 
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Fig. 2. — Observation from 2010 August 15. Top: nuclear spec- 
trum (region A of Figure [TJ. Middle: stellar population obtained 
from averaging the spectra from regions B and C in Figure [T] Bot- 
tom: nuclear spectrum after the subtraction of the stellar popula- 
tion. 

evidence for a significant contribution of an AGN contin- 
uum, as can be seen in Figure [51 The effect of a nuclear 
continuum, if present, would be to dilute the stellar ab- 
sorption lines, making them shallower at the nucleus than 
outside. The subtraction of the stellar population tem- 
plate, as described above, after normalization to the flux 
in the continuum of the nuclear spectrum would result 
in the presence of "fake" emission features in the resid- 
ual spectrum and a possible slope (in the case of a blue 
continuum, for example), which are not seen. 

We have performed a sequence of tests in order to esti- 
mate the contribution from a featureless continuum (FC) 
that could be "hidden" in the nuclear spectrum as fol- 
lows. We added the contribution of the FC, assumed in 
the form Fi, oc i/^-^, to the stellar population template 
described above, at increasing percent contributions of 
its flux at 5800A. The resulting composite (FC-|-stellar 
population) spectrum was then renormalized to the same 
flux at 5800A, and we subtracted the stellar population 
template from it. We concluded that for the short wave- 
length range and the signal-to-noise ratio of our obser- 
vations, in order to have clear residuals showing fake 
emission features and a blue continuum we need the FC 
continuum to contribute with at least 15%-20% of the 
flux at 5800A. We thus can say that a possible FC con- 
tribution to the nuclear spectrum is lower than 20% of 
the continuum flux at 5800A. We consider this value an 
upper limit, as a comparison between the stellar contin- 
uum from our spectra and the AGN continuum obtained 
in the work of Nemmen et al. (2006) shows that the stel- 
lar continuum in our aperture is about 30 times stronger 
than the AGN continuum in the optical. 

3. RESULTS 

After subtracting the contribution of the stellar popu- 
lation, the flnal step before the analysis of the double- 
peaked proflle was to normalize the spectra using as 
reference a previous flux-calibrated spectrum. We 
used as refere nce the spectru m of 1991 November 2 
(|Storchi-Bergmann et al.l [l993[ ). and assumed that the 
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flux of the narrow lines Ha, [Nii]AA6548, 6584, and 
[Sil] AA6717, 6731 did not vary during the last 20 years. 
This is a common assumption in this type of study, as the 
narrow lines originate in extended regions with spatial 
scales much larger (~ 100 pc) than that of the accretion 
disk (~ 10^'^pc). Thus, their response time to variations 
in the ionizing flux is correspondingly much larger than 
that for the line-emitting accretion disk, and cannot re- 
spond to variations of the ionizing flux that are as fast 
as those we have detected. 

The normalized spectra from our 2010-11 observa- 
tions are compared with the reference spectrum of 1991 
November in FigureO (pairs of spectra are shifted from 
each other in the vertical direction for clarity). In this 
figure we show the three spectra from 2011 January, the 
three from 2011 February, and the two from 2011 March, 
respectively, on top of each other because we did not find 
any significant variation within each group. 

Note that the double-peaked Ha profile from 2010 Au- 
gust is very different from that of 2010 March. In 2010 
March, the red peak was stronger than the blue peak, 
while in 2010 August, the blue peak was stronger than 
the red. This inversion of the relative intensity of the 
peaks occurred on a timescale of 5 months. In addition, 
in 2010 December, the profile showed a marked increase 
in its integrated flux compared to that of 2010 August, 
while the red and blue peaks showed similar strengths. 

From 2010 December 22 to 2011 January 11 we ob- 
served a '^ 25% decrease in the integrated flux of the 
double-peaked proflle. The three observations of 2011 
January were obtained within a time interval of seven 
days, and showed no significant variations. 

In 2011 February, the three observations were obtained 
within a time interval of six days, and again we did not 
see significant variations among the three profiles. How- 
ever, the total flux continued to decrease compared with 
that of January, with the blue and red peaks showing 
similar strengths. 

In 2011 March, we obtained two observations separated 
by a time interval of four days, between which we again 
did not find significant variation in the double-peaked 
profile. But the total flux decreased further, and the 
profile showed the red peak somewhat stronger than the 
blue peak. 

In summary, this series of observations obtained be- 
tween 2010 March and 2011 March reveal: (1) a marked 
increase of the broad-line flux between 2010 August and 
2010 December; (2) a significant decrease of the flux on 
timescales longer than a week; (3) no significant changes 
in the profile on time intervals shorter than a week; (4) 
an inversion in the relative strength of the blue and red 
peaks on a timescale of 5 months. 

3.1. Measurements of the double-peaked profile 
characteristics 

We have quantified the observed variations by mea- 
suring the wavelengths of the blue and red peaks of the 
profiles, their peak fiux densities, and the integrated flux 
of the broad line, as described below. We list the results 

in Table El 

In order to isolate the broad double-peaked Ha 
emission and perform the measurements, we flrst 
evaluated the contribution of the narrow lines 
Hanarrow A6564, [N ii] AA6549, 6584, [Sii] AA6718,6732 



and [O i] AA6302, 6365 by fitting gaussians to these lines 
and then subtracting them. The flux of the double- 
peaked line (Fbroad) was measured by integrating the 
total fiux under the profile after this subtraction. The 
average uncertainty in i^broad is 1 x 10""'^* erg cm~^ s~^. 

The wavelengths of the blue and red peaks (As and Xji) 
and the corresponding peak flux densities (Fb and Ffi) 
were obtained by fitting three Gaussians to the double- 
peaked profile, one for the blue peak, another for the 
red peak and a third Gaussian for the center of the pro- 
file. Fb and Fr are the peak flux densities of the Gaus- 
sians flttcd to the blue and red peaks, respectively, while 
Xb and A/j are the corresponding wavelengths. We have 
tried other methods as well, such as using splines to fit 
the peaks, but concluded that the fit of Gaussians gave 
equivalent but more robust results. The average uncer- 
tainties in the measurements are « 4Afor Xb and Ai^ 
and sa 6 X 10"-^^ erg cm~^ s~^A~^ for Fb and F^. 

From Ab and Xr we obtained the blue and red peak ve- 
locities Vb and Vr relative to the systemic velocity of the 
galaxy, adopted as the one corresponding to the narrow 
component of the Ha emission, as we have done in our 
previous studies (SB03). Figure 0] shows the temporal 
evolution of Vb and Vr. The velocities of the blue peak 
range between —5500 and —3500 km s~^ while for the 
red peak they are between 4000 and 6000 km s^^ . A com- 
parison between these peak velocities with those from the 
period 1991-2001 (SB03) shows that, in the new observa- 
tions, the blue and red peaks always have larger peak ve- 
locities than the minimum values observed between 1991 
and 2001, which are « —2500 km s~^ for the blue and 
« 3250 km s~^ for the red, as indicated by the lower gray 
horizontal bars in both panels of FigurelH Most velocities 
are instead similar or higher than the maximum values 
observed between 1991 and 2001, w —4500 km s^^ for 
the blue and ~ 4500 km s~^ for the red peak. These are 
indicated by the upper gray horizontal bars in FigurelH 
The values of Vs and V^ were obtained by using the 
expression for the relativistic Doppler effect. 

Figure [5] shows the variation of the integrated broad 
Ha emission line fiux from 2010 March to 2011 March. 
The integrated flux is always above the minimum value 
observed in the period 1991-2003 (lower gray bar in the 
figure), being close to the maximum value from 2010 
March to 2011 January (upper gray bar in the figure). 
A maximum in the fiux occurred between 2010 August 
15 and 2010 December 22, when the integrated flux rose 
above the maximum value observed between 1991 and 
2001. After 2010 December 23, there was an abrupt de- 
crease in the broad-line flux, by w 25% of its maximum 
value in 20 days. On the other hand, it can be observed 
that there is little variation on timescales shorter than 
a week (the maximum time interval between the obser- 
vations obtained in the same month). Thus signiflcant 
variations in the broad-line flux occur on a timescale just 
above a week. From Figure EJ one can derive a decrease 
of w 12% of the flux in 10 days. 

Figure [6] shows the temporal evolution of the peak 
fluxes of the blue and red peaks as well as that of their 
ratio Fb/Fr. On 2010 March, the red peak was stronger 
than the blue, while in 2010 August the blue peak was 
stronger. Thus, within a time interval of 5 months, there 
has been an inversion of the Fb/Fr ratio, and, in the 
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following 7 months this ratio decreased to a value which 
is just below 1, while in 2010 March this value was 0.7. 
Thus the change from a stronger red peak to a stronger 
blue peak was faster than the change back to a stronger 
red peak. From Figure[6] we thus derive a timescale for 
the inversion of the relative intensity of the two peaks of 
«5 months. 

Figure [7] shows that there is an inverse correlation be- 
tween the velocity separation of the red and blue peaks 
Vr - Vb and Fbroad, observed both in our recent obser- 
vations, represented by the circles, as well as in the pre- 
vious observations of SB03, represented by the crosses. 
This inverse correlation between Fbroad and Vr— Vb can 
be interpreted as due to: (1) when the central ionizing 
source becomes brighter, the line emission from the outer 
radii become relatively more important, as the radiation 
of the ionizing source reaches larger distances in the disk, 
where the velocity is lower; (2) when the central source 
is fainter, the inner radii are favored, as the ionizing ra- 
diation reaches smaller radii, where the velocity is larger. 

A comparison of the two data sets in Figure[7] shows 
that our recent data define a sequence parallel to that 
of the previous data, displaced from the previous data 
due to larger velocity separations between the two peaks 
and to higher fluxes. The larger velocity separation of 
the two peaks suggests that the line emitting portion of 
the disk is closer to the ionizing source or has a larger 
contribution from gas closer to the source. 

4. ACCRETION DISK MODELS 

In order to investigate the origin of the observed varia- 
tions in the double-peaked profile, we have used accretion 
disk models to reproduce the profile. Figure |6] shows that 
the double-peaked profile still shows variations in the rel- 
ative intensities of the two peaks, alternating between a 
stronger blue peak and a stronger red peak, as previously 
observed by SB03. This pattern of variation supports the 
origin of the profile in an accretion disk whose emission 
is not axisymmetric. In SB03, the relative strengths of 
the two peaks and their variation were better explained 
by emission from a circular accretion disk with a single 
spiral arm. The enhanced emission in the arm, as it ro- 
tates in the disk, produces the observed asymmetry of 
the profile. We will thus also adopt in our present mod- 
eling the single spiral arm model to fit the double-peaked 
profiles. 

The source of radiation which ionizes the disk in or- 
der to drive the Ha emission is also relevant to the 
modeling. As pointed out in the Introduction, in the 
case of NGC 1097, this source is an RIAF internal t o 
the line-emitting part of the disk ("N emmen et al .1120061 ). 
[Dumont & CoUin-Souffrin ( 1990) have found that, in the 
case of an accretion disk illuminated by a non-thermal 
continuum source with a spectral energy distribution ex- 
tending up to the hard X-ray band, there is a "radius 
of maximum emission" or "saturation radius" , at which 
the Ha emission reaches a maximum value. Thus we 
consider an additional model in which the disk emissiv- 
ity may increase with radius up to the saturation radius 
and then decrease outward. 

With the above considerations in mind, we then fit- 
ted the double-peaked Ha line profiles with two different 
models, described below. In both models we assume a 
circular emitting disk, or wide ring (as it has a hole in 
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Fig. 3. — Nuclear spectra from 2010 March to 2011 March each 
compared with the reference spectrum of 1991 November 2 (gray). 
For 2011 January-March we show the spectra obtained a few days 
apart plotted on top of each other, as no significant changes were 
found between them. 

the middle), with a perturbation in the emissivity in the 
form a spiral arm. In the first model which we call Spi- 
ral Arm Model the variations in the relative strengths 
of the peaks are controlled by the rotation of the spiral 
arm while the separation of the peaks is controlled by an 
emissivity law which is a simple power-law function of the 
radius. By changing the power-law index we can weight 
the emissivity of the disk toward inner regions (leading 
to a broader profile) or the outer regions (leading to a 
narrower profile). In the second model which we call the 
Saturated Spiral Model we have changed the emissivity 
law to a broken power-law, which can produce a ring of 
maximum emissivity between the inner and outer radii. 
In this case, the separation of the peaks is controlled by 
the radius of this ring of maximum emissivity. 

4.1. Spiral Arm Model 



We adopt the formulation described in [Gilbert et al.l 
(|1999| ). in which: (1) the line emitting portion of the disk 
is circular and located between an inner radius ^i and an 
outer radius ^2 (where ^ is the disk radius in units of the 
gravitational radius rg = GM,/c^, c is the light speed, G 
is the gravitational constant and M, is the mass of the 
black hole); (2) the disk has an inclination i relative to 
the line of sight (zero degrees is face on). Superimposed 
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TABLE 2 
Measurements of the Profile Properties 



Date 



(A) 



(A) 



Fb"-^ 



(10 



-15j,j.g j,j^ 2 



^A^i) (lO"^"' erg cnr 



2 



iA-1) (10- 



-'^ broad 
-15 ^j,g c]-Q~^ 



-^) 



2010 Mar 4 


6483. Oil. 6 


6656.5±0.5 


0.760±0.06 


2010 Aug 15 


6469.0±1.0 


6671. 7±2. 2 


0.760±0.06 


2010 Dec 22 


6474.5±2.5 


6659.6±2.6 


1.160±0.02 


2011 Jan 11 


6462. 1±2. 8 


6667.0±1.9 


0.830±0.01 


2011 Jan 15 


6463.8±3.2 


6667.7±1.9 


0.840±0.01 


2011 Jan 18 


6463.7±1.4 


6669.3±1.2 


0.860±0.02 


2011 Feb 8 


6453.5±2.1 


6677.3±1.2 


0.550±0.03 


2011 Feb 10 


6456.7±2.3 


6675.3±2.3 


0.600±0.05 


2011 Feb 14 


6460.1±2.2 


6675.4±2.2 


0.570±0.01 


2011 Mar 17 


6463.8±1.4 


6678.2±1.0 


0.490±0.01 


2011 Mar 21 


6463.8±0.7 


6678.9±1.1 


0.470±0.01 



0.970±0.02 
0.570±0.03 
1.070±0.02 
0.700±0.02 
0.710±0.01 
0.730±0.01 
0.580±0.02 
0.600±0.03 
0.580±0.02 
0.520±0.03 
0.500±0.02 



197.4±7.7 
186.1±4.4 
258.5±7.2 
210.1±6.2 
216.7±4.0 
218.7±1.2 
153.2±1.5 
162.7±1.8 
150.0±5.2 
129.2±3.0 
127.1±3.2 



'^ \b and \ji are the peak wavelength of the blue and red peaks respectively. 

^ Fg and .F^ are the corresponding peak fluxes. 

'^ Abroad is the integrated flux of the broad double-peaked Ha emission line. 
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Fig. 4. — Evolution of the peak velocities of the blue (top panel) 
and red (bottom panel) peaks from 2010 March to 2011 March. 
The gray horizontal bars show the minimum and maximum values 
(plus uncertainties) in the spectra from 1991 to 2001 (SB03). 

on the axisymmetric emissivity of the circular disk, there 
is a perturbation in the form of a spiral arm. We adopt 
only one spiral arm, as we have already verified that two 
or more arms do not reproduce well the variations in the 
profile (SB03). 
The total emissivity of the disk is given by: 



e(e,0)=e(C)^l + |cxp 



4 In 2 



S^ 



■-^of 



A 
--exp 



41n2 



(27r - (/) + iPof 



where 



eiO = eor 



(1) 



(2) 
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Fig. 5. — Evolution of the integrated flux of the broad double- 
peaked Ha line Fbroad from 2010 March to 2011 March. The gray 
horizontal bars show the minimum and maximum values (plus un- 
certainties) in the spectra from 1991 to 2001 (SB03). 



is the axisymmetric emissivity of the disk, A is the bright- 
ness contrast between the spiral arm and the underlying 
disk, and the expression between square brackets repre- 
sents the decay of the emissivity of the arm as a function 
of the azimuthal distance 4> — ipo from the ridge line to 
both sides of the arm, assumed to be a Gaussian function 
with FWHM S (azimuthal width). 

The relation between the azimuthal angle 0o a-nd the 
angular position ijjq of the ridge of emissivity on the spiral 
arm is given by 



^po = (l)Q + 



tanp 



(3) 



where (f>o is the azimuthal angle of the spiral pattern, p 
is the pitch angle and ^gp is the innermost radius of the 
spiral arm. 

The specific intensity from each location in the disk, 
in the frame of the emitting particle is calculated as: 



I{^j(P,^e) 



477 (27r)l/2cr 



(4) 
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Fig. 6. — Evolution of the peak flux densities of the blue peak 
(top panel) , the red peak (middle panel) and of the ratio between 
the two (bottom panel) from 2010 March to 2011 March. 
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Fig. 7. — The inverse correlation between Fbroad and Vr — Vb- 
Circles represent the measures for this work and crosses represent 
the data from SB03. The inverse correlation for the new observa- 
tions however is displaced to larger velocities and fluxes relative 
to previous observations (SB03) . This suggests that emission from 
the inner parts of the disk became more important relative to that 
of the outer parts. 
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Fig. 8. — Best flts of the spiral arm model for the epochs from 
2010 March 4 to 2011 March 21, obtained by varying (/iq (orienta- 
tion of the spiral pattern), A (contrast between the spiral arm and 
the disk), and q (the slope of the emissivity law). The values of 
the parameters are listed in Table |3] 

where ly^ is the emission frequency, vq (Ha 6564. 6A) is 
the re st frequency and a is th e local "broadening param- 
eter" (IChen fc HalperniriQSl) . 

In the fits of the individual profiles we first used the 
same set of parameters of SB03, namely: the inner and 
outer radii of the line emitting portion of the disk ^i — 
^sp — 450 and ^2 = 1600, the inclination angle i — 34°, 
broadening parameter a = 1200 km s~^, pitch angle p = 
50° and S = 70°. In order to reproduce the evolution of 
the profile we varied the parameters A, q and (Jjq- The 
variation of (j)o regulates the evolution of the shape of the 
profile, while a variation in A, together with 00, regulates 
the relative intensities of the two peaks. The value of 
q regulates the radius in the disk of the maximum line 
luminosity, which corresponds to the inner regions for the 
broadest profiles, or to the outer regions for the narrowest 
profiles. 

The best-fitting models are compared to the observed 
line profiles in Figure HI while the corresponding param- 
eter values {A, q and (j>o) are listed in Table |3l The 
contrast of the spiral arm A ranged from 2.0 at the first 
epoch, decreased to a minimum value of 0.5 in 2011 Jan- 
uary, and increased again to 2.0 in 2011 March, while q 
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TABLE 3 

Parameters for the Spiral Arm Model 

Fits. 









(Ao 


Shift 


Date 


A 


9 


(deg) 


(km s-l) 


2010 Mar 4 


2.0 


-2.0 


130 





2010 Aug 15 


1.0 


1.0 


250 





2010 Dec 22 


0.3 


1.0 


330 





2011 Jan 11 


0.5 


1.5 


340 


-460 


2011 Jan 15 


0.5 


1.5 


340 


-460 


2011 Jan 18 


0.5 


1.5 


340 


-460 


2011 Feb 8 


1.0 


2.0 


360 


-460 


2011 Feb 10 


1.0 


2.0 


360 


-460 


2011 Feb 14 


1.0 


2.0 


360 


-460 


2011 Mar 17 


2.0 


2.0 


400 


-460 


2011 Mar 21 


2.0 


2.0 


400 


-460 



increased from the first to the last epochs. An abrupt 
variation in q occurred between 20f0 March and 2010 
August when the value of q changed from —2.0, favor- 
ing the outer parts of the disk, to 1.0, favoring the inner 
parts of the disk. The (/>o values were constrained to 
provide a monotonic rotation of the spiral pattern, from 
130° to 400°. The spiral arm has thus almost completed 
one revolution in one year. 

For the spectra from 2011 January to 2011 March, an 
improvement of the fit can be obtained by allowing a 
small, ad hoc blueshift of the central wavelength of the 
profile relative to the wavelength of the narrow Ha line 
(6564.6A). This blueshift is w -460 kms-\ as hsted in 
the last column of Table [31 The improved fits with the 
blueshifts are shown as dashed lines in Figure [5] 

In Figure [HI we show the evolution of the (pQ values ob- 
tained from the fits, represented by the circles in the fig- 
ure, while the black solid line is the best linear regression 
to the points and the gray band covers the uncertainty 
in the regression. The linear regression has angular co- 
efficient (j)o = 0.680 ±?0. 02 day~^, which represents the 
angular velocity of the arm and implies a rotation period 
of F « 17.7 ±0.5 months. 

4.2. Saturated Spiral Model 

In an effort to improve the fit to the broad double- 
peaked Ha profiles, we modified the spiral arm model 
to allow a more flexible radial emissivity law. We 
have chosen an emissivity law which allows the exis- 
tence of a "radius of maximum emissivity" between the 
inner and outer radii of the disk at which the broad 
Ha line reaches a maxim u m int ensity, as proposed by 
iDumont fc: Collin- SouffrinI (|1990( ). The emissivity in- 
creases with radius until the saturation radius and then 
decreases outward. 

The radius of maximum emissivity was incorporated 
by changing the axisymmetric part of the emissivity of 
the disk (Equation [2]) as follows: 



e(0 



eo6 



-(«l-92)£- 



, Cl < C < Cg 
, Cg < C < 6 



(5) 



where: the new parameter ^q is the radius of maximum 
emissivity, or saturation radius, at which the emissivity 
law changes; gi is the index of the emissivity law for the 
regions internal to C? (Ci < C < Cg); I2 is the index for 
the regions external to Cg (Cg < C < C2); Ci and ^2 are 
the inner and outer radii, respectively. 
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Fig. 9. — Evolution of the (/>o values obtained from the fit of the 
spiral arm model. The solid line is a linear fit to the data and the 
shaded bar covers the uncertainty in the fit. 



TABLE 4 

Parametes for the Saturated Spiral 

Model Fits 



Date 



<j)0 Shift 

(deg) (km s^-"-) 



2010 Mar 4 


3.0 


1200 


140 





2010 Aug 15 


1.0 


800 


220 





2010 Dec 22 


0.2 


1000 


330 





2011 Jan 11 


1.0 


650 


340 


-500 


2011 Jan 15 


1.0 


650 


340 


-500 


2011 Jan 18 


1.0 


650 


340 


-500 


2011 Feb 8 


1.0 


550 


360 


-550 


2011 Feb 10 


1.0 


550 


360 


-550 


2011 Feb 14 


1.0 


550 


360 


-550 


2011 Mar 17 


1.5 


600 


400 


-230 


2011 Mar 21 


1.5 


600 


400 


-230 



We kept the value (72 = 3.0, as discussed by 
IDumont fc Collin-SouffrinI ()1990f ). who showed that be- 
yond the saturation radius the flux emitted by the line 
decreases as oc S^~^ . After many tests, gi = —2.0 gave 
the best fits, and was fixed at this value. The negative 
value allows the increase of the emissivity until the sat- 
uration radius. The parameter S,q was allowed to vary. 
We kept the parameters of the spiral arm at the values 
Csp = ?ij P = 50°, 5 = 70°, as well as the inclination of 
the disk at i = 34°. 

We have varied the broadening parameter and found 
that the best value to reproduce the shape of most pro- 
files was a — 900 km s~^. 

In summary, in the Saturated Spiral model, the only 
parameters we needed to vary in order to reproduce the 
evolution of the profile were ^q, A, and 0o- The changes 
in the separation between the blue and red peaks were 
obtained by changing ^q] when ^^ is smaller, the separa- 
tion between the blue and red peaks increases. 

The best-fitting models are compared with the ob- 
served line profiles in Figure llOi and the corresponding 
parameters values (^g, A and </)o) are listed in Table S) 
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Fig. 10. — Best fits for the hybrid model, obtained by varying 
<j>o (orientation of the spiral pattern), A (contrast between the disk 
and the spiral arm) and ^q (radius of maximum emissivity). The 
values of the parameters are listed in Table [4] 



The parameter ^q ranged from a maximum value of 1200 
in 2010 March when the separation between the peaks 
was the smallest, to a minimum value of 550-600 in 2011 
February-March, when the separation between the peaks 
was the largest. The parameter A assumed the value 
3.0 in 2010 March, when the red peak was significantly 
stronger than the blue, as seen in FigurefTOl while in 2010 
December this parameter reached the minimum value of 
0.2, when the peaks had similar strengths. Similarly to 
the previous model, 4>o varied from 140° in the first epoch 
to 400° in the last. 

As in the case of the Spiral Arm model, from 2011 Jan- 
uary to 2011 March, the fit was improved by introducing 
a blueshift which is listed in the last column of Table |31 
This blueshift ranges from « —230 to —550 kms~^, and 
the improved fits with the blueshifts are shown as dashed 
lines in Figure [TOl 

In Figure [11] we show emissivity maps constructed us- 
ing the parameters of the best fits using the Saturated 
Spiral model. In 2010 March the radius of maximum 
emissivity ^g was closer to ^2 , when the outer parts where 
brighter, while in 2011 March ^g was closer to ^1 and the 
inner parts were brighter. 



5. DISCUSSION 

5.1. Timescales 

In order to discuss the physical processes driving the 
changes in the profiles, we first review four relevant 
timescales for the accretion disk in NGC 1097: the 
viscous timescale (tvjsc), the sound-crossing timescale 
(ts), the dynamical timescale (rdyn) and the light 
travel tim,es c ale (t i). These timescale are presented in 
iFrank et al.l ()2002[ ) and can be modified in order to ex- 
press the d isk radius in terms of the gravitational radius 
()Eracleou£|[l99l iLewis et al.llMoh : 



io°m:'^£,. 



3/2^5/4 -4/5 



a_ 



'M 



-3/10 

_i years 



70Ms^3T, 



1/2 
5 years 



'^dyn 



SMsCs^^^ months 



n = 6Ms^3 days 



(6) 



(7) 
(8) 



(9) 



where ^3 = lO""^^, Mg = M. x IQ-^Mq^', a = lO^^ x 
a-i is the Shakura-Sunyaev viscosity parameter, T^, is 
the temperature in units of 10^ K, and M_i is the mass 
accretion rate in units of 0.1 Mq yr~^. 

For NGC 109 7, M, = 1.2 x IO^Mq 

(jLewis fc EracleousI 120061) . and M/Mwm ^ 0.0064 
(in Eddington units. iNemmen et al.l (J2006D ). corre- 
sponding to M = 0.018 M0yr~^. Adopting typical 
values of a = 0.1 and T5 = 0.1, a characteristic radius 
for the disk is ^ = 1000. Adopting the above values of 
M, and M, we obtain the following estimates for the 
timescales: 

1- T'visc ^ 10^ years: this large value for the vis- 
cous timescale does not support that the observed 
changes in the profile — on timescales from weeks to 
years — are due to inward transfer of mass due to 
viscous torques; 

2. Ts ~ 265 years: this is also much larger than the 
observed variations timescales; 

3. Tdyn <^ 7 months: the dynamical timescale depends 
on the radius, with the inner parts of the disk hav- 
ing shorter timescales than the outer parts. Adopt- 
ing ^3 — 0.45 for the inner radius, as obtained 
from our modeling, the corresponding dynamical 
timescale is just above 2 months. The new ob- 
servations presented here show that the variation 
timescale of the relative strengths of the blue and 
red peaks {Fb/Fb) is w 5 months (see Figure [6]), 
thus compatible with the dynamical time. In our 
modeling, this variation was reproduced by the ro- 
tation of the spiral pattern around the disk, whose 
period is close to the dynamical time at ^2; 

4. T; > 7 days: considering the range of radii of the 
disk, this timescale ranges from «4 to wl2 days. 
This is the shortest timescale of the disk. Ac- 
cording to our data (see Figure [5]), the smallest 
timescale of the profile variations — namely varia- 
tions in the integrated fiux of the double-peaked 
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line — occur in <7 days. This supports the inter- 
pretation that the variability in the integrated flux 
-Fbroad occurs in the hght travel timescale, as a re- 
sult of variations in the ionizing flux. Similarly, 
the variations in Vb and Vr — see Figure 0] — which 
correlate with Fbroad, also occur in the light travel 
timescale. According to the models, these are due 
to changes in the region of maximum emissivity of 
the disk. 

5.2. Constraints on the physical mechanisms causing 
the variations 

The timescales discussed above provide constraints on 
the mechanisms that cause the changes in the double- 
peaked profile. Our data clearly reveal two distinct 
timescales: (1) the first corresponding to changes in the 
integrated flux of the double-peaked profile, which is just 
over 7 days; and (2) the second corresponding to the 
changes in relative strengths of the blue and red peaks, 
which is in the range 5-6 months. 

We have identified the first timescale with the light 
travel timescale. In order to model the corresponding 
variations in the flux of the double-peaked line, and at 
the same time reproduce the variation of the width of 
the profiles it was necessary to vary the disk emissivity. 
When the flux is higher, the profile is narrower, and the 
emissivity is higher in the outer parts of the disk, while 
when the flux is lower, the proflle is broader and the 
emissivity is relatively higher in the inner parts of the 
disk. In the Saturated Spiral model, what varies is the 
radius of maximum emissivity. 

The physical scenario we favor is the following: the 
double-peaked profile originates in gas rotating in a 
thin a ccretion disk, which is ionized by an RIAF 
(jNemm en et al. 2006) located inside the inner radius of 
the line-emitting portion of the disk. As the variations 
occur on the light travel timescale, our interpretation is 
that the variation in Fbroad can be considered a "rever- 
beration" effect produced by variations in the ionizing 
flux from the central RIAF, and allows a test for the val- 
ues of the radii obtained via the models. Calculating the 
values for the inner and outer radii of the disk in light 
days we obtain for ^i « 3 light days and for ^2 ~ H light 
days, which is consistent with the variability timescale of 
« 7 days that we observe. 

In order to further confirm the reverberation hypothe- 
sis, it would be necessary to do a proper "rever beration 
mapping campaign" (e.g., IPeterson et al.l (J2004D ) of the 
accretion disk via simultaneous observations of the RIAF 
emission (e.g., in X-rays) and the double-peaked profile. 
We plan to pursue such a monitoring campaign in the 
future, including a more extended spectral coverage to 
include the blue part of the optical spectrum which will 
also allow us to test other models for t he origin of th e 
continuum, such as the one proposed bv iGaskelH ()20lH ). 

We have identified the second timescale — of the varia- 
tions in the relative strengths of the blue and red peaks — 
with the dynamical timescale. This means that the asym- 
metry in the profile is due to a non-axisymmetric emit- 
ting structure - modeled as a single spiral arm - rotating 
in the disk, with a rotation period of ~ 18 months. This 
timescale and its interpretation differs from our previ- 
ous one (SB03). In that paper we concluded that the 
variation of the relative strengths of the two peaks oc- 



curred in the sound-crossing timescale. This distinct in- 
terpretation was due to the fact that the period we have 
derived then for the variation was 5.5 yr, and at that 
time our estimate for the SMBH mass was w 10^ Mq, 
giving a smaller value for the sound-crossing timescale, 
consistent with the derived period of variation. With 
a revised SMBH mass two orders of magnitude larger, 
the revised sound-crossing timescale is much larger than 
the timescale of the variation of the relative fluxes of the 
blue and red peaks we have now derived from our more 
frequent observations. 

Is important to point out that the perturbation does 
not need to be a spiral arm, even though we used the 
spiral arm prescription to model it. In the images of 
the disk, constructed using the best model parameters 
for the different epochs, shown in Figure [TTl the spiral 
arm actually looks more like an azimuthal arc. There- 
fore, it could be that this arc is actually a clump that 
was prod uce by self-gravity and the n sheared (see dis- 
cussion in lFlohic &: EracleousI ()2008[ )). Another alterna- 
tive is that this feature is an irrad iation-induced warp 
(|Pringle|[T996t IMalonev et"aIlll996D . But although the 
irradiation-induced warp model can r eproduce the ob - 
served variability of the line profiles (|Wu et al.l l2008[ ) , 
the required rotation timescale is incompatible with the 
mass of the SMBH in NGC 1 097 (see Equation (4) in 
iStorchi-Bergmann et al.l ()1997D V On the other hand, the 
sheared clump would orbit the SMBH on the dynami- 
cal timescale and therefore could reproduce the 5 month 
variation that we observe. 

Figure[7] shows that the profiles got broader and 
brighter in the recent (2010-2011) observations when 
compared to the previous ones (1991-2001), although at 
both epochs there is an inverse correlation between the 
flux and width of the profiles. This inverse correlation is 
expected in the reverberation scenario discussed above, 
but there is an offset between the previous and present 
correlations. A comparison between the parameters of 
the present and previous spiral arm model shows that the 
emissivity law in the recent epochs favors a higher contri- 
bution of the inner regions relative to the outer regions 
of the disk, leading to a broader profile. One possibil- 
ity to explain this increase in emissivity of the innermost 
regions is an increase in the gas density in these regions 
of the disk in the most recent epochs when compared to 
the period 1991-2001. 

A comparison between Figures l8l and ITOl shows that the 
Saturated Spiral model gives a better fit to the double- 
peaked profiles. In order to quantify the goodness of 
the fits, we have obtained the root mean square (rms) 
deviation between the observed profiles and the models. 
In almost all observations the Saturated Spiral model 
provided a smaller rms when compared with that for the 
Single Spiral arm model, confirming what is seen in the 
fits. 

Finally, we point out that, although the models can 
in general reproduce the Ha profiles, there are several 
epochs in which we observe an excess flux in the blue 
wings of the profiles relative to the models. This ex- 
cess flux is observed from 2010 August to February 2011 
(Figures 151 and ITUl) . In addition, from 2011 January to 
March we observe a systematic blueshift of the center of 
the line, from —230 to —550 km s~^ which could be due 
to an accretion disk wind. Therefore, it will be worth- 
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while in future work to compare the observed line pro- 
files with mod e ls tha t include an outflowing wind (e.g., 
IHall fc Chaiitl !^0W )) 

5.3. The mass of the supermassive black hole 

The physical value of the disk timescales, such as the 
light-travel timescale (r/), as well as the inner and outer 
radii of the accretion disk, are strongly dependent on the 
mass of the SMBH. Thus, uncertainties in its determi- 
nation must be taken into acc ount. We have a dopted 
here the SMBH mass value of L ewis fc Eracleoui (2006): 
M, = 1.2 ±0.2 X 10^ Mq. 

The range of values for the inner and outer radii of 
the accretion disk (^i — 450 and ^2 = 1600, respec- 
tively) resultin g from the u ncertainties in the SMBH 
mass quoted bv lLewis fc Erac lcous (2006) are 2.7 < ^1 < 
3.6 light-days and 9.1 < ^2 < 12.8 light-days. But 
iLewis fc EracleousI (|2006|)~obtain the value for the SMBH 
mass using the M, — a relation of lTremaine et al.l ()2002t ). 
pointing out that they do not take into account the scat- 
ter in this relation. They consider only the uncertainty in 
the coefhcients of the relation that results from the scat- 
ter (effectively they assume that the observed scatter is 
not intrinsic). If we now consider that the uncertainty 
in the SMBH mass is also affected by the sc atter in the 
M, — a relation from iTremaine et al.l ()2002[ ). which can 
introduce an uncertainty in the mass of ~ 0.3 dex, the 
resulting minimum and maximum values for the inner 
and outer radii become: (1) for M, w 0.6 x 10® Mq, 
the minimum radii are ^1 ~ 1.5 light-days and ^2 ~ 5.5 
light-days; (2) for M, « 2.4 x 10* Mq, the maximum 
radii are ^1 ~ 6 light-days and ^2 ~ 22 light-days. 

Thus, even considering the maximum uncertainties 
which could affect the determination of the SMBH mass, 
the resulting values for the inner and outer radii of the 
disk are still of the order of light-days. The above re- 
sults support the interpretation that the shortest varia- 
tion timescale of 7 days is the light-travel time between 
the ionizing source at the center of the disk and the mean 
radius of the disk. In fact, this observed timescale can 
actually constrain the values of f 1 to < 7 light-days and 
of ^2 to > 7 light-days. 

In addition, if we assume that 7 light-days is the mean 
distance between the broad-line-emitting clouds and the 
SMBH, we can estimate the mass of the SMBH inde- 
pendently from the model, as follows. If this system of 
clouds is flattened, as suggested by their double-peaked 
profile, and we adopt as the corresponding Keplerian ve- 
locity around the SMBH the mean velocity of the peaks — 
« 4650kms~i— we obtain Af, = 1.1 x 10* Mq, for an 
inclination of 34° relative to the plane of the sky (as 
obtained from the disk model). Th is value is in good 
agreem ent with the one obtained bv ILewis fc EracleousI 
(|2006l) . 

6. CONCLUSIONS 

We report the discovery of short timescale variations 
in the broad double-peaked Ha emission line from the 
LINER nucleus of NGC 1097, observed in 9 epochs be- 
tween 2010 March and 2011 March. They comprise: (1) 
variations in the integrated flux and width of the line over 
a timescale of > 7 days, which we identified as the light- 
travel time between the ionizing source and the disk; it 
is the first time that such short timescale variations have 



been seen in a double-peaked profile from the nucleus of 
a galaxy; (2) variations in the relative intensity of the 
blue and red peaks on a timescale of 5-6 months, which 
is compatible with the disk dynamical timescale. 

Using two distinct accretion disk models we reproduced 
these variations through the combination of two effects: 

1. short timescale changes in the disk emissivity — 
namely in the contrast A and index q in the Spiral 
Arm model and A and S.q in the Saturated Spiral 
model; 

2. the rotation of the spiral arm — via the change of 
the orientation 0o of the arm in both models — 
resulting in a estimate for the rotation period of 
the arm of « 18 months. 

The values of the inner and outer radii, obtained from 
previous modeling, were kept in the present modeling and 
are, respectively, ^1 = 450 and ^2 = 1600 (corresponding 
to w 3 and 11 light-days, respectively). These values are 
in good agreement with the short timescale variations 
seen in our observations. This agreement supports the 
identification of the shortest timescale variations with 
the light travel time between the ionizing source and the 
accretion disk, in one hand, and gives further support 
to the model and the derived parameters, on the other 
hand. 

Our recent observations support the scenario in which 
an accretion disk whose emissivity changes due to rever- 
beration of a variable ionizing continuum (presumably 
from an RIAF), has a non-axisymmetric feature rotating 
in the dynamical timescale. 
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Fig. 11. — Images showing the disk emissivity from 2010 March to 2011 March, using the hybrid model. White represents the brightest 
regions, and the observer is to the bottom. The disk parameters for these epochs are listed in Table l4l 
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